fdatool
G
SOS

Hd=filter1;
x=randn(1,1000);
y=filter(Hd, x);
plot([x; y]')

sx=20*log10(abs(fft(x)));
sy=20*log10(abs(fft(y)));
figure(2);
plot([sx; sy]') 

